Emergence of pulled fronts in fermionic microscopic particle models 
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We study the emergence and dynamics of pulled fronts described by the Fisher-Kolmogorov- 
Petrovsky-Piscounov (FKPP) equation in the microscopic reaction-diffusion process A + A ^ A 
on the lattice when only a particle is allowed per site. To this end we identify the parameter 
that controls the strength of internal fluctuations in this model, namely, the number of particles 
per correlated volume. When internal fluctuations are suppressed, we explictly see the matching 
between the deterministic FKPP description and the microscopic particle model. 
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The study of diffusion-limited reaction processes has 
shown the important role of internal or microscopic fluc- 
tuations in low dimensions Q, Q 1^ - Mean- field approx- 
imation for those processes assume that diffusive mix- 
ing is much stronger than the influence of correlations 
produced by reactions. However, diffusive mixing is not 
strong enough in low-dimensional systems and fluctua- 
tions might modify the dynamics or induce nonequilib- 
rium phase transitions. While this behavior is observed 
in different situations, there is an special interest in the 
problem of front propagation in reaction-diffusion sys- 

tems iiiiaiilBlIllIlllllIllIllIllia In 

this paper we concentrate in microscopic lattice reaction- 
diffusion models whose mean- field approximation is given 
by the FKPP equation ilj| 
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where p{x, t) is the local concentration of particles. Such 
an equation display traveling-wave solutions of the form 
p = p[S^) with ^ = x — vt which invade the unstable phase 
p ~ from the stable phase p = ^1/^2 £ind travel with 
velocity v > vq — 2y/Dki. For steep enough initial con- 
ditions, the solution selected for large times is the one 
with minimal velocity, Uq, which is known to be a pulled 
front, since it is essentially "pulled along" by the growth 
and spreading of small perturbations in the leading edge 
where p <C 1 |3j 0- Microscopic fluctuations are ex- 
pected to modify macroscopic properties of pulled fronts 
at two levels: (i) first because the deterministic descrip- 
tion breaks down at small densities p ^ 1/N where N 
is the number of particles, which introduces an effective 
cutoff in the FKPP equation. Due to the importance of 
the tail development in pulled fronts, several front fea- 
tures are dramatically affected by this effective cutoff [sl : 
for example, the selected velocity converges as 
to the mean- field value vq, (ii) second, because internal 
fluctuations are present and could interplay with or even 
destroy pulled front development and dynamics [slF^ITsl. 
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One of the most studied microscopic models is the 
reversible reaction model A A + A on a lattice 
[ill 11 Ills 111 111 El 111 111 111 113. Inthebosonic 
version of this model |3 , the number of possible particles 
per site is unbounded and thus the balance between birth 
and coagulation gives an average number of particles per 
site N. If iV ^ 00 the reaction is well stirred within each 
site and the front dynamics is described by the mean- 
field approximation For very large N discreteness 
effects remain and produce the predicted velocity correc- 
tion V — Vq ^ ln~^ N 11^3- In the case of the fermionic 
version of the A ^ A + A model only a particle is allowed 
per site. The main reason to consider exclusion is that, 
for some values of the parameters, the model is analyti- 
cally tractable HHHIlJI and/or simulations are easier 
than in the bosonic version. Exact results are available 
for the two interesting regimes in the model: reaction- 
limited regime 14] ■ where coarse-grained density front 
profiles are described by the mean-field FKPP equation, 
and diffusion-limited regime H H Q in which internal 
fluctuations dominate front propagation and the mean- 
field approximation is not valid. Our purpose in this 
paper is to put these two results in a general framework 
that can describe the emergence of pulled fronts in this 
fermionic model. This is done by identifying the control 
parameter that modulates the effect of internal fluctua- 
tions on the front propagation model. As we will see, this 
parameter also controls the development of the tail front 
and establishes the appearance of pulled fronts. 

In the A <-> A+A model in one dimension, particles are 
allowed to occupy lattice sites and can undergo the fol- 
lowing moves: (i) Diffusion to any one of its two neighbor 
lattice sites with a diffusion rate D (ii) Birth: occupied 
sites spontaneously generate particles at neighbor lattice 
sites with rate p, (ii) Coagulation: a particle can get an- 
nihilated with death rate rj if one of its two neighboring 
filled lattice sites is occupied. The fermionic nature of 
the model makes diffusion and birth only possible if the 
neighboring lattice site is empty. The mean-field descrip- 
tion of this model is given by the FKPP equation with 
ki = 2p, k2 = 2(p + r]). Starting from an initial condition 
in which occupation number is only different fron zero on 
the right side of a site, a front develops and advances as a 
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FIG. 1: Velocity as a function of the birth rate fi for different 
values of the coagulation rate rj. Simulations are done with 
D = 0.5. Dashed lines correspond to the predictions v — fi 
0,13 and V = vo 14|. Inset: Velocity corrections as a function 
of /i for the case = 0. Dashed line is the power law N'^^'"^ 
while dash-dotted line is the prediction 7r^/(21n^ A''*) of 0. 



FIG. 2: Schematic diagram of condition ((21 in the rj, n phase 
space. Shaded area correspond to those values of the param- 
eters for which A'^* > 1. Solid lines correspond to the solution 
of Eq. Q for different values of A''*. Dashed lines are the 
following sets of parameters: ?) = -D[^I3:A' = ^ [Ti| and 
rj = 0.002 used in Fig.Cl 



function on time. Operationally, the front position Xf{t) 
is determined by a local average of density of particles 
over intervals of length = 2D/vq, which is the width 
of the deterministic front selected for our initial condition 
by the FKPP equation 19]. Thus, front position is 
the point where this coarse-grained density equals Pq/2 
[T^ fisll . Other definitions of the front position yield the 
same results 0| . Since results only depend on the ratios 
D/jj, and D/rj we set D = 1/2 throughout this paper. 
After a transient time (which could be long), the front 
advances linearly, i.e. {xf{t)) — vt, where (• • • ) stands 
for average over different realizations. In this regime, 
statistical properties of the front with respect to the nor- 
malized coordinate ^ = x — vt, are independent of time. 

In this model, there are two special cases for which 
exact results are available: when rj = D the model is 
solvable usiiig the method of inter-particle distribution 
functions P, 'a] . In that case fronts advance with velocity 

V — fi which shows how internal fluctuations can domi- 
nate the system behavior in this diffusion-limited regime. 
On the other hand, when 77 = it was proved in [l3 | 
that fronts approach asymptotically the FKPP equation 
predictions {v = vq) in the limit D/ fj, 00 (reaction- 
limited regime) , while v = D + fi in the opposite regime 
D/fj, ^ (diffusion-limited regime). In Fig. ^ we show 
the results of our simulations for the velocity of the front 

V as a function of /i for different values of ri. Our re- 
sults are consistent with the exact results 0, IE Q and 
previous simulations of this model • For an interme- 
diate case < T] < D we observe that, for some values 
of n, the velocity seems to approach the deterministic 
value Vq. However, for small enough value of /x internal 
fluctuations seem to dominate and the velocity deviates 
strongly from vq. 

To understand this behavior, let us recall that mean 



field approximation in the A A + A is only valid 
when diffusive mixing is strong enough. Specifically, this 
happens when the typical distance traveled diffusively by 
a particle between reaction events, Ijj, is much larger that 
the typical distance between particles, Ip 3\. In that case, 
particles are well stirred within cells of size Id and thus, 
mean-field approximation is valid for the coarse-grained 
density of particles over cells of size I d as shown in |Q| . 
In our model we have that Id = min(A/ Df/j,, yjD/rj), 
while we approximate Ip by the average distance between 
particles in the stable phase Ip = Pq^ = {n^-rf)/ p.. Since 
we are interested in the propagation of pulled fronts, 
which are only driven by the birth term, our condition 
to approach the mean-field approximation is then given 
only by /i: 



(2) 



Interestingly, this condition is equivalent to N* = 
A~^po ^ 1) where N* is approximately the number of 
particles within an interval of length A~^. Thus condi- 
tion (0) also means that the number of particles within 
the typical length scale of the front (its width) is large. 
Our results in Fig. ^ are then easily explained in terms of 
condition Q : when N* ^ 1 internal fluctuations should 
be unimportant within cells of size Id — and front 
propagation should approach asymptotically the FKPP 
predictions. Since 



N* = 



(3) 



then A^* > 1 only happens when {Dli]f D / fi :$> I for 
fixed values of D and rj. We show the condition N* > 1 
in Fig. 121 along with the different set of parameters used 



in Fig. n and in other works 0, 0, 0| • Outside the re- 
gion N* > 1 internal fluctuations dominate and fronts 
are not described by the FKPP equation. This is the 
case and r] = fj, ^Jl- the intermediate 

case < rj < D we can have values of /i for which N* 
is relatively large and fronts seem to approach the deter- 
ministic value of vq , which explains the behavior observed 
in Fig. ^for rj = 0.002. Note however, that although be- 
ing in the region TV* > 1 is the minimum requirement for 
our model to approach the mean-field description ^ a 
finite value of A'^* means that fronts are still subject to 
internal fluctuations and discreteness effects which pro- 
duce a (strong) correction to the velocity. Only in the 
limit N* — > oo do both effects become negligible and the 
A ^ A + A system is effectively described by the FKPP 
equation. This is the case for ry = 0, — > [20|. 

An interesting question is whether N* plays any role 
like the average number of particles per site, N, in 
bosonic models JfeiSiildl- In those models, it is observed 
that the deterministic description of a pulled front given 
by the FKPP equation is valid until the density drops to 
p ~ N^^ which produces an effective cutoff in the tail 
of the front and modifles its velocity To check this 
possibility, we have measured in our simulations the av- 
erage distance of the last particle from the front position, 
which is observed to saturate to a constant value for 
long enough times. It is obvious that for ^ > ^* the 
continuum description of the front breaks down and we 
expect this to happen when there is only a particle in each 
coarse-grained site of length A^-'^, i.e. when p{£,*) ~ aX, 
where a is a constant. When internal fluctuations are 
irrelevant, i.e. when N* 3> 1, we assume that the con- 
tinuum description is still valid up to f * and taking 
that p ~ pctX£,e^^^ for > 1 for a pulled front 2^ we 
obtain 
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(4) 



Solutions of this equation for ^* with N* given by Q are 
compared with our simulations in Fig. O We see that 
for A^* > 1 Eq. 10} gives a rather accurate prediction 
of This corroborates our assumption that a pulled 
front described by the FKPP equation develops even for 
moderate values of N* up to the point where p ~ {N*)~^. 

An important consequence of Eq. Q is that N* con- 
trols not only the size of internal fluctuations but also 
the appearance and length of the tail in the pulled front. 
Thus, when iV* < 1 internal fluctuations dominate and 
also the tail length is roughly zero, A^* ~ 0. This means 
that the front is basically a shock wave with height pQ. 
Actually, the exact solution when = -D shows that 
the particles behind the leading one remain distributed 
as in the stable phase (p = po) at all times, which con- 
firm our picture (see Fig. ^ . This shock wave shape is 
also observed in the case rj ^ D when p ^ rj, where it 
is found that A^ ~ 0. In the intermediate case in which 
77 ^ we see that the front develops a tail which is de- 
scribed by the FKPP equation only for a given interval 
of values of p (see Fig. |31for = 0.002). 




FIG. 3: Position of the last particle from the front position, 
^* , as a function of p from simulations. Dashed lines are Eq. 
@ with a = 1.4. Inset: occupation probability for the site 
just before the last particle (circles). Solid line is given by 
Eq. with V taken from simulations and dashed line is the 
approximation p* ~ a/A used to get Eq. |0J 2^ . 



In the case 77 = we have studied the correction to 
the velocity as a function of N* and observe in the inset 
Fig.^that it decays like {vq —v)/vo ~ (iV*)~^/^ which is 
consistent with simulations of other microscopic bosonic 
models H, 01 for moderate values of number of particles 
N. This results stresses the equivalence of N* with the 
role that the number of particles plays in other micro- 
scopic models. In particular we expect the correction to 
be {v — vo)/vo ln~^ N* for very large values of N* [2ll| . 

Since the last particle is, on average, at a certain dis- 
tance from the front position^ their velocities coincide. 
This fact was used in [Til [T5| to estimate the velocity 
of the front by counting possible forward and backward 
hopping rates: 



v = p~p*ir]-D) 



(5) 



where p* is the probability of having a particle behind 
the last one. Several approximations can be made for 
the value of p* [TlllTsl . For example, in is was taken 
as p* ~ Po, i.e. p* is given by the probability to find 
a particle in the stable phase. Clearly, this approxima- 
tion is only valid in the case in which fronts are like a 
shock-wave, i.e. when N* ~ 1 because then the last par- 
ticle is very close to the stable phase. In the case in 
which a pulled front develops {N* ^ 1) we find that Eq. 
(O still holds: since the last particle is on average at a 
fixed distance from the front position, we can approxi- 
mate p* ^ pqXS,*^ '^^ = a/A which is the concentration 
of particles at 1^*. Our simulations for 77 = confirm the 
validity of this approximation (sec Fig. |3Jl which brings 
out the effective matching between the continuum de- 
scription given by the FKPP equation for ^ < ^* and the 
microscopic character of the model for C C*- 

Another interesting property of front propagation is 
the wandering of the position of the front around its mean 
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FIG. 4: Time evolution of the diffusive spreading of the front 
position for different values of ^ (solid lines) and r) = 0: fj, — 
10"^ 10"^, 10"*, lO"'^ from top to bottom. Inset: Diffusion 
coefficient as a function of A'^* for r] — 0. Dashed line is the 
power law (iV*)"^''^. 



value A2(t) = {{xf{t) - {xf{t)))'^). Several studies of dif- 
ferent microscopic models shows that A^(t) — 2D ft for 
long times and that the diffusion coefficient Df 

depends on the number of particles N. Specifically, it 
was found that Df ^ N^^^^ for moderate values of N 
dllll, while Z?/ ~ ln~^ iV for very large values of iV |3|. 
Our simulations for the A + A A system shows that 
fronts move diffusively for all values of the parameters. 
In the case 77 = 0, in which the model approaches asymp- 
totically the FKPP equation, we get Df ~ {N*)~^/^ (see 
Fig. like in bosonic models ^3 which stress once again 
the fact that N* plays the role of the number of par- 
ticles in this fermionic model. Moreover, we found for 
small times that as N* increases, the correlation between 
the time development of the front and internal fluctua- 



tions produces superdiffusive motion of the front position 
A2(<) ~ t'^" with ly ~ 0.8. Once the front tail is devel- 
oped (which happens at t ~ k^^), front position starts to 
wander diffusively. Finally, our results for the diffusion of 
the front indicate that as the front approaches the deter- 
ministic FKPP equation, internal fluctuations make the 
front move diffusively at times t > k^^, independently of 
N* . We do not observed any signs of subdiffusive behav- 
ior conjectured by some authors |4, 12] for pulled fronts 
subject to noise. This supports the idea that pulled fronts 
subject to internal noise belong to a different universality 
class than those subject to external noise p^. 

In summary, we have identified the parameter that con- 
trols the strength of microscopic fluctuations for the front 
propagation problem in the fermionic model A A + A, 
namely the number of particles N* per coarse-grained 
site of length A^^. When A^* ^ 1, internal fluctuations 
are suppressed and the front becomes a pulled front like 
those of the FKPP equation Moreover, our stud- 
ies about the length of the tail, the velocity of the front 
and its diffusion show that A''* plays the same role as the 
number of particles in other microscopic bosonic models. 
Finally, it is interesting to note that in the A ^ A + A 
model, the velocity of a macroscopic object such as the 
front is related to the microscopic motion of the last par- 
ticle, something also observed in other works 8, IX, 1^ . 
We hope our results will help to understand the dynam- 
ics of fronts in microscopic fermionic reaction-diffusion 
models and its relevance when discussing properties of 
the FKPP equation subject to internal noise |3.ll2.ri3|. 
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